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The Landau-gauge gluon propagator is computed for the SU(3) gauge theory on lattices up to 
a size of 32 3 x 200. We use the standard Wilson action at (3 = 6.0 and compare our results with 
previous computations using large asymmetric and symmetric lattices. In particular, we focus on 
the impact of the lattice geometry and momentum cuts to achieve compatibility between data from 
^ | symmetric and asymmetric lattices for a large range of momenta. 
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I. INTRODUCTION AND MOTIVATION 



The investigation of the Landau-gauge gluon propagator 



in lattice QCD dates back to more than 20 years. Despite the efforts made by a number of authors, there are 
questions which remain to be answered. For large momenta, let us say q > 2 GeV, the results from different groups 
provide a consistent picture of the propagator and agree well with calculations performed using other non-perturbative 
1—1 1 techniques. On the other hand, the infrared limit is still an open issue and has been a field of intense research in the 
last years (see the contributions to this conference, e.g. and references therein). From the point of view of lattice 
simulations, the questions to be answered yet are numerous and sometimes are even not easy to handle with. For 
' example, it is still under debate how to deal with the Gribov ambiguity in lattice simulations, and how close present 
lattice data come to Landau-gauge gluodynamics in the continuum and infinite- volume limit. 

In order to access the infrared limit of the gluon propagator, two of us (O.O. and P.J.S.) recently proposed and 
explored the use of large asymmetric lattices @, H, i-e- L 3 x T with T ^> L. The price of relying on such kind 
of lattices are the control, or the lack of it, of additional finite-volume effects coming from a breaking of the Z4 
t-H ' symmetry, a remnant of the 0(4) continuum symmetry on a symmetric hypercubic Jattice. When previous studies on 
C**"* . symmetric lattices have shown strong finite- volume effects in the infrared region [f|, the situation is more dramatic 
for asymmetric lattices. For example, there the gluon dressing function Z(q 2 ) = q 2 D(q 2 ) computed at equal time-like 
^ , and spatial momenta are not necessarily compatible within pure statistical errors in the low-momentum region. 

On the other hand, the access to very low momenta is much more computationally intensive in simulations on 
symmetric lattices compared to those on asymmetric ones. Therefore, if somehow the asymmetry-induced finite- 
volume effects were brought under control, data at much lower momenta than currently available could be obtained. 
Indeed, then still the infinite volume limit has to be taken, but in a situation where more data were available in the 
infrared momentum region. 

Having now access to a considerably larger spatial volume for the asymmetric case, in this study we report on some 
first results obtained comparing data on symmetric and asymmetric lattices, namely 16 3 x 256, 18 3 x 256, 32 3 x 200 
and 32 . In particular, we look for regions in the lattice momentum space where the differences between time- like 
and spatial momenta disappear and where not. For our simulations we use the standard Wilson gauge action with 
j3 = 6.0 fixed. This value corresponds to an inverse lattice spacing of about a -1 = 1.94 GeV. To relate our data at 
the different lattice momenta to their continuum counterparts we use 

2 / nn \ 

qtl = a \l) ' rfy = 0, 1, 1, (2) 



where is the lattice extent in direction fi. Definitions and details on the gauge fixing are given in [7| and for the 
32 3 x 200 data in [8(. In the following, whenever possible, a Z3 average over equivalent momenta is performed. 
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Bare Gluon Dressing Function 

32 3 x200, |i = 6.0 




Bare Gluon Dressing Function 

32 J x200. p = 6.0 




FIG. 1: The gluon dressing function on a 32 x 200 lattice. The different momentum cuts are discussed in the text. 



II. ASYMMETRIC LATTICE: 32 3 x 200 

The gluon propagator and dressing function were computed for 39 configurations using a 32 3 x 200 lattice. The 
data shows discretization effects similar to those seen in the symmetric lattices. Indeed, it is well-known that for 
momenta above q > 1 GeV the propagator is not a simple function of q 2 = q 2 ^ alone. The traditional approach is to 
apply momenta cuts which reduce the dependence of the propagator and the dressing function on other Z\ (in our 
case Z-i) invariants to a unique curve. This is better seen in the gluon dressing function. In order to illustrate this 
effect, in Fig. Q] we plot the gluon dressing function for different choices of momenta (purely time-like and different 
cuts of spatial momenta). 

The plot shows that, within our limited statistics, there is very good agreement between the dressing functions 
computed using purely spatial on-axis and purely temporal momenta. The figure does not include the dressing 
function for all purely spatial momenta. However, in what concerns the purely spatial momenta, the gluon dressing 
function for on-axis momenta evolves typically along the lower edge of the spatial (including off-axis) momentum 
data. The diagonal choice of momenta, i.e. the cylindrical cut Q where n M w ±n v (see the left plot in Fig.[T]) picks 
up an unique propagator which is slightly above the propagator for the on-axis choice for momenta at q larger than 
~ 1 GeV. Note that this "democratic" choice of momenta has been successfully used to suppress discretization effects 
such that data from different volumes and lattice spacings match at larger momenta. 

In the right plot of Fig. [1] we also show two naive generalizations to the asymmetric case of the on-axis case and of a 
"democratic" choice of momenta (the cylinder cut), labelled as (111X) and (1116) cuts. The former includes momenta 
only of type n M = (±1, ±1, ±1, ±n t ) with n t = 0,1, ... , T/2 (softening of on-axis momenta), while the latter includes 
momenta defined around the direction n^ = (in, in, in, ±6n t ) with n = 0, 1, — 1 and n t = n, n ± 1, n ± 2. 

This direction is close to the diagonal in the elongated volume (remember that T/L = 6.25). Note that the (111V) 
cut reproduces the results of the symmetric lattice for the on-axis choice of momenta, while the (1116) cut follows the 
symmetric lattice data for the cylindrical and conical cuts. Given that the cylindrical and conical cuts seem to reduce 
the finite- volume effects for momenta above 1 GeV, Fig. [1] suggests that for the asymmetric lattices one should use 
the (1116) cut, or a variant of it, for that momentum region. 



III. SYMMETRIC LATTICE: 32 4 



The gluon dressing function Z(q 2 ) = q 2 D(q 2 ), computed for an ensemble of 50 gauge configurations on a 32 4 lattice, 
is given in the left plot of Fig. [5] applying various cuts. The lattice data shows similar discretization effects as for the 
asymmetric lattice. This is illustrated in the right plot of the same figure. There, data for the two lattices 32 4 and 
32 3 x 200 are shown for two different momentum cuts at larger momenta and good agreement is found. Moreover, 
within our limited statistics, the cuts produce similar effects for both lattices where the on-axis data lie in both cases 
systematically below cylinder-cut data for q > 2 GeV. 

In Fig. [3] the lattice gluon dressing function for a 32 3 x 200 and a 32 4 lattice are compared for momenta below 2 
GeV. Again, two different cuts (on-axis and cylindrical) are considered. We find that in this momentum range the 
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Bare Gluon Dressing Function 

32 4 vs. 32 3 x200, p = 6.0 
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FIG. 2: The gluon dressing function for a 32 4 lattice is shown on the left hand side. The right figure compares the same data 
with those on a 32 3 x 200 lattice. Different cuts have been applied to the data. 

Bare Gluon Dressing Function 
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FIG. 3: The gluon dressing function at low momenta for the 32 and the 32 x 200 lattices for two different momenta cuts. 



dressing functions Z{q 2 ) for the two lattice geometries are in good agreement, even though we cannot compare at the 
low-lying momenta. A comparison with data from larger symmetric lattices is necessary to become more confident in 
this. At least, the good matching between 0.5 GeV and 2 GeV is encouraging in what concerns the use of asymmetric 
lattices to extract reliable infrared properties in future lattice simulations. 



IV. THE IMPACT OF THE SPATIAL VOLUME 



Now we discuss the volume dependence of the propagator D(q ^ 0) as function of the physical momentum q. 
In Fig. [4] we compare the gluon propagator on one hand and the dressing function on the other hand for various 
asymmetric lattices, for 32 4 (all at (3 = 6.0) and for the continuum Dyson-Schwinger solution of Ref. [10]- All 
propagators were renormalized according to the condition 

£(* 2 W = ^ (3) 

with the choice \i = 3 GeV. For asymmetric lattices, the data is for time-like (on-axis) momenta. For the symmetric 
lattice, the plot includes only on-axis momenta. Fig. 0] shows that the two results become closer as the lattice volume 
increases. 
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FIG. 4: The gluon propagator (left) and the gluon dressing function (right) for various asymmetric lattices and for 32 4 at 
/3 = 6.0 . compared with the DSE solution. 
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0.5320t3o 
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TABLE I: Parameters corresponding to a fit of the lattice gluon dressing function according to Zi(q 2 ). q m ax and A are given 
in MeV. 



V. FITTING THE IR GLUON DRESSING FUNCTION 

In previous investigations it was verified that the lattice gluon dressing function and the continuum Dyson-Schwinger 
solution are well described not by a pure power law but by both functions 

ftaO ■ W) - "wrUvr (4) 

for momenta below ~ 700 MeV. 

The results of fitting the lattice dressing function for purely temporal momenta with Zi and Zu are reported in 
tables U and HH respectively; q max is the highest momentum included in the fits. Note that, the exponent k agrees 
within one standard deviation for the two largest volumes. Furthermore, for these lattices, n agrees with the estimate 
of O.O. and P.J.S. k ~ 0.53, from using ratios of propagators to suppress the volume dependence. The ratio 
method discussed in [3j, if applied to the 32 3 x 200 data, estimates n — 0.565 ± 0.040 for the infrared exponent. 

The results give k consistently above 0.5. If this really represents the infrared asymptotics, it supports a vanishing 
q — > limit of the gluon propagator D(q =/= 0). Moreover, one should keep in mind that the fits to a pure power law 
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TABLE II: Parameters corresponding to a fit of the lattice gluon dressing function according to Zu(q 2 ). q m ax and A are given 
in MeV. 
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provide always a k < 0.5, with k increasing with the lattice volume. It should be noted, however, that so far lattice 
simulations have always reported a finite and not vanishing gluon propagator at zero momentum [12j . 

VI. RESULTS AND CONCLUSIONS 

In this work the gluon propagator and dressing function has been analysed for various asymmetric lattices and a 
comparison to 32 4 data has been done. Despite the observed finite- volume effects, for volumes as large as 32 3 x 200 the 
dressing function Z(q 2 ) for purely temporal momenta agrees well, within the available statistics, with the corresponding 
function for purely spatial on-axis momenta. 

In what concerns the 32 3 x 200 and 32 4 data, the momentum cuts produce similar results for the full range of 
momenta, i.e. the gluon propagator/dressing function for on-axis momenta are systematically below cone-cut or 
cylinder-cut data for q > 1 GeV. 

The behaviour of the lattice infrared gluon dressing function Z(q 2 ) is well described by the two ansatze Zj and Zjj 
for q < 700Af eV. The fits to the data provide n values which support a vanishing zero momentum limit of the gluon 
propagator for all the lattices reported here. The measurement of the infrared exponent k for the two larger lattices 
suggests a value k ~ 0.53, in agreement with the estimate discussed in jH. 
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